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Abstract 

The spectroscopic factors for the low-lying quasi-hole states observed in the 
16 0(e, e'p) 15 N reaction are reinvestigated with a variational Monte Carlo cal- 
culation for the structure of the initial and final nucleus. A computational 
error in a previous report is rectified. It is shown that a proper treatment of 
center-of-mass motion does not lead to a reduction of the spectroscopic fac- 
tor for p-shell quasi-hole states, but rather to a 7% enhancement. This is in 
agreement with analytical results obtained in the harmonic oscillator model. 
The center-of-mass effect worsens the discrepancy between present theoretical 
models and the experimentally observed single-particle strength. We discuss 
the present status of this problem, including some other mechanisms that may 
be relevant in this respect. 
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I. INTRODUCTION 



In present-day (e, e'p) experiments [lj one can determine quite accurately single-particle 
(s.p.) overlap functions between an A-particle nuclear target in its ground state and low- 
lying (quasi-hole) states of the (A — 1) residual system. The spectroscopic factor, defined 
as the norm of the overlap function, is equal to unity only in the fixed-center mean-field 
approximation. Experimentally one finds that spectroscopic factors of about S ~ 0.6 — 0.7 
are needed to explain the data. Such deviations of S from unity are normally ascribed to 
nucleon-nucleon (NN) correlations. The observed values indicate the sensitivity of the quasi- 
elastic (e, e'p) cross section to correlation effects in the initial and final state. Several nuclear 
many-body calculations have been made in attempts to explain the reduced spectroscopic 
factors of the pi/ 2 and P3/2 quasi-hole states seen in the 16 0(e, e'p) reaction 

Fixed-center mean-field wave functions are not realistic for light systems like 16 O, since 
they are not translationally invariant, i.e. they contain spurious center-of-mass (cm.) mo- 
tion. Variational Monte Carlo (VMC) calculations of the 16 ground state ]7| and the 
15 N(p 3 / 2 ) quasi-hole state , in which the nuclear wave functions are explicitly transla- 
tionally invariant, have been made to address this problem. It was reported in || that cm. 
correlations reduced the p-shell spectroscopic factor by ~ 12% even in absence of dynam- 
ical NN correlations. However, as was noted in ||, this reduction contradicts the results 
obtained in the harmonic oscillator (h.o.) model, which predicts an enhancement by 7% 
for the p-shell spectroscopic factor in 16 O |§ due to cm. correlations. The discrepancy 
with the h.o. model values indicated an error in the computer program of H]. It has now 
been identified, and the corrected results are presented in the third Section of this paper. 
The correction simply involves a rescaling of the s.p. overlap function reported in [[3] by 
the factor (y§) 3//2 . This brings the results without dynamical NN correlations into perfect 
agreement with the h.o. results, but the full calculation now predicts a p-shell spectroscopic 
factor S p = 0.98, considerably different from the experimental values. 

The remainder of the paper is organized as follows. In Section II we examine the concept 
of overlap functions with translational invariance of the nuclear wave functions taken into 
account, and point out some consequences for the description of knock-out reactions. We 
also mention general results of the h.o. model and apply it to the case of the 16 nucleus. The 
theoretical description of the quasi-hole states observed in the 16 0(e, e'p) reaction is treated 
in Section III, where we correct the results of 0. It is pointed out that a proper treatment 
of cm. motion in other calculations would also worsen the discrepancy with experimental 
data. In Section IV we look at the order of magnitude of two simple corrections to the 
(e, e'p) cross section. The present status of the problem is discussed in Section V. 

II. CENTER-OF-MASS MOTION AND KNOCK-OUT REACTIONS 

A. Overlap functions in self-bound systems 

For systems which are localized around a fixed force center (e.g. the electrons of an atom), 
one defines a s.p. overlap function between normalized A and (A — 1) particle systems as 

4>{x A ) = VA / dxi . ■ ■dx A „ 1 $\ A _ 1) {x 1 , . . . ,x j4 _ 1 )$( j4) (x 1 , . . .,x A ). (1) 
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(The notation Xj includes the spatial coordinate, r^, and the appropriate spin and isospin de- 
grees of freedom.) In the case of Fermi systems, the <&(a) an d $(a-i) are antisymmetric wave 
functions, and as a consequence, the normalization of the overlap function (or spectroscopic 
factor S) has the property 

S = j ' dx A \<P(x A )\ 2 <1. (2) 

In nuclei there is no fixed external force center, but the nucleons are localized around 
their cm. due to their mutual interactions. The eigenstates of such a self-bound system can 
be factored as 



®(A)(xi, ...,xa) = exp(iK • Ka)^(a){xi, ■ ■ -,xa), 



(3) 



with ^(A) the intrinsic, translationally invariant, wave function of the system, and the plane- 
wave factor describing the cm. motion; K is the total momentum and is the position 
of the cm. The properties of the overlap functions are somewhat different in this case. A 
detailed analysis of many s.p. quantities in self-bound systems (the one-body density matrix, 
natural orbitals, spectral function etc.) will be given in a future publication [ |lC)|j . 

For the present discussion it will be sufficient to note that, if the intrinsic wave functions 
are normalized according to 



J dxi . . . dx A S(R A )\^ {A)(xi, ■ ■ . , x A )\ 2 = 1, 



then the overlap function tp must be defined as [11 



1p(x A ) = VA J dXt... dxA-lS(R.A-l)^\A-l)( x li ■ ■ • ,xa-i)^(A)(xi 



,x A 



(4) 



(5) 



As shown in Section II. B, the overlap function defined according to Eqs.(g-jg) is indeed the 
counterpart for self-bound systems of Eq. ([[]), and it is the natural quantity appearing in the 
description of knock-out reactions when some standard approximations are made. 

It should also be kept in mind that the spectroscopic factors in self-bound Fermi systems, 



S 



dx A \ip(x A )\' 



(6) 



can be larger than one [with deviations of the order 0(1/ A)}. Although this was already 
known before PJT2||, the fact does not seem to be widely appreciated. The most extreme- 
though somewhat forced-example is that of two identical fermions (e.g. two spin-up neu- 
trons). The intrinsic wave function can be written as ^(2) = — r^), with / dr\f(r)\ 2 = 1. 
Then the overlap function with the (only) intrinsic one-particle state, \&(i)(ri) = 1, is given 
by ip(r 2 ) = y/2f(— r 2 ), and the corresponding spectroscopic factor is S — 2, regardless of 
the Fermi or Bose nature of the particles. A more realistic illustration for three particles is 
considered in the Appendix. 
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B. Intrinsic transition matrix elements for knock-out reactions 



In this Section we extend the standard derivation of the plane-wave impulse approx- 
imation (PWIA) amplitude in fixed-center systems to the case of knock-out reactions on 
self-bound systems, in order to justify the definition Eq.fl5|) for the overlap function. We 
consider for simplicity the transition matrix elements of a scalar one-particle external probe, 
e.g. the operator p(q) = X/^ =1 exp(iq- r^), between translationally invariant initial and final 
states of the form (|3|). The matrix element is readily separated into a momentum conserving 
5-function and the intrinsic transition matrix element, 

($ (A)/ |p(q)|$ (A)i ) = (2vr) 3 ( 5(K l + q-K / )(vI/ (A)/ |p(q)|vI/ (A)j ), (7) 

given by 

A 

(*(A)/|p(q)|*(A)i) = / cfei . . . dx A 8(R A )w ( ^ A)f (x 1 , . . . ,x A )^2exp(iq ■ r^^A^xi, . . .,x A ). 

(8) 

In knock-out reactions where the residual (A — 1) nucleus is left in a bound state ^m-i), 
the intrinsic final state can be specified by its asymptotic behavior, 

lim ^u^ p (xi,... J x A ) = Cexp[ip-(r A -K A -i)]^(A-i)(xi J ...,XA-i), (9) 

|rA-R.A-l|->°° 

with p the relative momentum between the knocked-out nucleon and the cm. of the remain- 
ing ones. 

In the plane-wave approximation, distortion effects on the escaping particle are neglected, 
and we simply approximate 

^(A)p(xi,---,x A ) « =^l{exp[ip- {r A -B. A -i)]^ff A -i)(x 1 ,...,x A -i)}, (10) 

^A\(A-1)\ 

with A the antisymmetrized sum over all coordinate permutations. The wave functions 
^(A) P are not orthogonal to the initial state ^( A )i, and also not mutually orthogonal, since 

(^ ( A)p|^(A)p') = J dxi . . .dx A 5(R A )^ A)p (xi, . . . ,x A )*(A) P '( a; i> ■■■,xa) (11) 
= 5(p-p')-C(p,p'), (12) 

with the non-orthogonality correction C(p, p') given by 

f 11 
C(p,pO = (27r)- 3 y dxdx' p(x,x') exp[ ir ■ (p' + ^—jp)]exp[-tr' ■ (p + ^—^p% (13) 

The quantity p(x,x') is the one-body density matrix of the intrinsic (A — l)-particle state 



p(x,x') = (A - 1) J dxi . . .dx A ^ 2 5(R A -2)^l A ^i)(xi, . . . , x A - 2 , x)*(a-i)(xi, • • .,x A - 2 ,x'). 

(14) 
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Similarly as for fixed-center systems, the nonorthogonality correction C(p, p') is suppressed 
by the momentum dependence of the Fourier transform of p(x, x') for momenta p or p' larger 
than typical momenta contained in the bound-state ^m-i), and can be neglected under the 
usual (e, e'p) kinematical conditions. 

The amplitude itself can now be written as 

A 

(# ( A) p |p(q) | %L)i) = I dx 1 ...dx A 5(K A ) V\ A _ 1) {x 1 ,...,x A -i) (15) 

j'=i 

exp[-zp • (r A - R A -i)] exp[zq • r,-] . . . ,£a), 

The impulse approximation implies that the nucleon which has absorbed the momentum of 
the probe gets ejected, and corresponds to retaining only the term with j = A in Eq . (|15|) . 
The result is 

r A — 1 
(*(A) P b(q)|*(A)i) ~ ApwiA(p,q) = {2tt)- 3/2 J dx A exp[-ir A ■ (p q))il>(x A ), (16) 

where ip(x A ) given by Eq.(|5|), i.e. the PWIA amplitude still scales with the Fourier transform 

of the overlap function, at missing momentum p j^q- 

The antisymmetrization correction terms (with j ^ A) in Eq.fJT^) can be recombined 
into 

-4A.s.(p,q) = (27r)~ 3/2 y dxdx'F(x, x') exp[ir • ( - q + - - - p)] exp[-zr ; ■ ( p + — )], 

(17) 

with 

F(x,x') = VA(A - 1) J dxi . . . dx A ^ 2 5CR A ^2)^l A ^i)(xi, . . . , xa-2, x)^( A )(xi, ■ • • , x A -2, x, x). 

(18) 

This correction is again suppressed by the momentum dependence of the Fourier transform 
of F(x, x') when p and q are large. The magnitude of the nonorthogonality and antisym- 
metrization corrections will be estimated in Section IV. B for the kinematical conditions used 



in the (e, e'p) experiment of fll"3|. 



C. Intrinsic Slater determinants 

Descriptions of all but the lightest nuclei usually involve fixed-center (shell-model) wave 
functions of the Slater determinant (SD) type, 

1 

$(A)(xi,...,x A ) = —=Det[(f) h:j (xi)] iij=li .._ tA . (19) 

These wave functions contain spurious cm. motion which does not correspond to the intrinsic 
degrees of freedom in self-bound systems. The correct intrinsic wave function describing the 
uncorrelated motion of A particles in A orbits around their cm. is given by 
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...,x A ) = -=Det[0 hj (ri - Ra, ai)kj=i,...,A, (20) 
v yi! 

where «j represents spin and isospin degrees of freedom. These intrinsic Slater determinants 
(ISD) are obviously translationally invariant, and should be used in combination with Eqs.(Q- 
d), to calculate s.p. overlap functions and spectroscopic factors. 

The effect of the spurious cm. motion on overlap functions and spectroscopic factors 
is non-trivial, even when considering transitions between a single ISD ^(a) like ( f20|) and a 
one-hole state ^>(A-\)h n obtained by removing the orbital h n , 

^{A-X)K ( x l> • • • ) x A-l) = I 1 Det[0 fe , (r, - Ra-1, Oii)]iJ=l,...,A;ijiA-Jjtn- (21) 

<{A-1)\ 



In contrast to the situation in fixed-center systems, the presence of the cm. 5-function in 
Eqs.(B-Ep makes the calculation of the normalization and overlap integrals quite complicated, 
even for ISD wave functions. Only if the orbitals that build up the ISD are chosen to be h.o. 
wave functions can the overlap function and spectroscopic factor be evaluated analytically. 
The simplest non-trivial illustration is for a h.o. well containing three neutrons, one spin- 
down and two spin-up; this is worked out in the Appendix. 

The general h.o. model has been studied in [0], where overlaps were considered between 
an A-nucleon ground-state configuration in the h.o. shell model, and the (A — l)-nucleon 
one-hole states. It was found that, in each spin-isospin space ar, the spectroscopic factor 
for the valence hole state (corresponding to a hole in the occupied shell with the largest 
oscillator quantum number N v ) is larger than one and given by 

/ A \ Nv 

^ T) = (^l) • ( 22 ) 

Moreover, for h.o. orbitals the spectroscopic sum rule remains satisfied in the intrinsic frame?]: 
in each spin-isospin space the sum of the spectroscopic factors of all one-hole ISD's 
yields the number of particles A^ T ^ 

En^Sj^Wt"), (23) 

h 

with n^ T ^ the number of particles occupying the main h.o. shell h. Since the sum rule is 
satisfied and the spectroscopic factor for the valence hole state is larger than one, the deeper 
hole states are partly spurious and have S < 1. 

For the case of 16 we have n s — 1, n p — 3, N v = 1; the h.o. model predicts for a hole 
in the p-shell S p — || ~ 1.07, and, through the sum rule (p^), S s = 4 — n p x i| = | for a 
hole in the s-shell. More realistic choices of the s.p. wave functions, such as Woods-Saxon 
(W.S.), preclude an analytical treatment of cm. effects. Nevertheless one can expect results 
close to the h.o. values for a light nucleus like 16 0, where h.o. and W.S. wave functions are 



In general, the spectroscopic sum rule is not satisfied by considering just one-hole states of ISD's 
made of non-h.o. orbitals; this is also illustrated for a three-neutron example in the Appendix. 
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rather similar. We checked this by a direct computation of the overlap functions in 16 O with 
W.S. wave functions in the ISD's, exploiting the fact that for ISD the many-body integrals 
in Eqs.(H-|^) can be reduced to a sequence of one-body integrals, e.g. we have (apart from 
normalization factors) 



A — 1 



x J rfkDet J cfxexp[zk • ^^.(^^.(r — — r A , a) 



• (24) 

i,j=l,...,A;i^n;j=im 



As expected the spectroscopic factors obtained with the W.S. wave functions practically 
coincide with the h.o. results. 

Note that the correlated many-body wave functions considered in the next section are 
more complicated than single ISD's, and Monte Carlo quadrature was used there to calculate 
the s.p. overlap functions and spectroscopic factors. We verified that, for the case of single 
ISD's, the Monte Carlo quadrature and a direct evaluation of Eq.(f24]) lead to the same 
result. 



III. SPECTROSCOPIC FACTORS FOR QUASI-HOLE STATES IN 16 Q 



In the variational Monte Carlo framework of IK?], the intrinsic wave function of the 16 
ground state has the form 



■ XA) 



F( Xl , 



.x A ) 



'A\ 



zDet 



i,j=l,...,A 



(25) 



The intrinsic Slater determinant is built up with the s.p. wave functions <ph, and incorporates 
the mean-field aspects. As explained in Section II. C, it is explicitly translationally invariant. 
The correlation factor F contains two- and three-particle correlations of central, spin, isospin, 
tensor and spin-orbit type, and is translationally invariant by itself. 

The correlation factor and the s.p. wave functions were determined by a minimization of 
the expectation value of the many-body Hamiltonian 



z ' nL i=i i<j i<j<k 



(26) 



where the Argonne v 14 model of the two-nucleon interaction |L4|| and the Urbana model VII 
of the three-nucleon interaction [|T3| were used. A detailed description of the variational 
procedure can be found in Q|. 

The intrinsic wave functions *&(A~i)hn °f ^ ne quasi-hole states in 15 N are approximated 

as 



ty( A -l)h n (Xl, ■ ■ -,xa-i) = F(x!, . . .,x A -i)- 



1 



1 A\ 



=Det 



i,j=l,...,A;i^A;j^n 



(27) 



i.e. by retaining the correlation factor and s.p. orbitals of 16 and omitting in the determinant 
the column with the corresponding <ph r 
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The overlap function for the p 3 / 2 quasi-hole state was then calculated according to Eqs. (fi- 
ll). When dynamical NN correlations were neglected (by putting the correlation factor 
F = 1), a spectroscopic factor S P3/2 = 0.88 was found. This case corresponds to the s.p. 
overlap between single intrinsic Slater determinants discussed in Section II. C, for which the 
harmonic oscillator model predicts S p = 1.07. The large discrepancy indicates an error in 
the computer program used in [Q], which has now been identified and corrected. 

To correct the error in 0], the s.p. overlap functions calculated with translationally 
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invariant wave functions simply have to be rescaled by (-jf) 2 - The corrected spectroscopic 



factor with only cm. correlations is thus given by S Pa/2 = 0.88 x (if) 3 = 1.07, and is now 
in perfect agreement with the h.o. model. The results including dynamical NN correlations 
in H must be changed likewise, and we get S P3/2 = 1.06 when central NN correlations are 
added, and S P3/2 = 0.98 for the complete calculation. The value S Ps/2 = 0.90, obtained by 
including dynamical NN correlations, but no cm. correlations, does not change. 

We note that the proper treatment of cm. motion makes the discrepancy with the 
experimentally obtained spectroscopic factors severely worse. In |H| a spectroscopic factor 
S Pl/2 (g.s.) = 0.61 is reported for the l/2~ ground state in 15 iV, whereas the lowest 3/2~ state 
at 6.32 MeV has ^^(6.32) = 0.53. The experimental low-lying p 3 / 2 strength is fragmented 
over three states at 6.32, 9.93 and 10.70 MeV, of which the 6.32 MeV state is the dominant 
one, with 87% of the total strength. Since we did not include fragmentation due to low- 
energy configuration mixing, our variational result S P3/2 = 0.98 should be compared to the 
total experimental value S P3/2 = 0.53/0.87 = 0.61. The cm. correction is not limited to 
variational calculations, but affects other theoretical models as well, such as the Green's 
function calculations in for which cm. motion was neglected. Thus, the result S P3/2 = 

0.91 in would (in first order) be changed by cm. correlations to S P3/2 = 0.91 x ^| = 0.97, 
a value similar to ours, while that of || (where low-energy fragmentation is taken into 
account) is increased from 5^3^(6.32) = 0.76 to S P3/2 (6.32) = 0.81. 



IV. OTHER CORRECTIONS TO THE (E, E'P) CROSS SECTION 

A correct treatment of cm. motion enhances the spectroscopic factor of valence hole 
states at the mean-field level, leading e.g. to a 7% enhancement for the p-shell spectroscopic 
factors in 16 O. This results in present calculations [|]-|| of the 15 N spectroscopic factors 
giving values that are much larger than the experimental values. In this Section we discuss 
two simple mechanisms, not taken into account in that could possibly further reduce 

the theoretical cross section for knock-out of a p-shell nucleon. Firstly, the spectroscopic 
factor itself could change at the mean-field level if we allow for different mean fields in the 
target and residual nucleus. Secondly, corrections to the PWIA reaction amplitude due to 
nonorthogonality and antisymmetrization may become more important when cm. motion 
is taken into account. 



A. Different s.p. orbitals in target and residual nucleus 

The assumption [implicit in Eq.(^7|)1, of having the same set of s.p. wave functions 
describing the target and residual nucleus, is likely to become less adequate as A decreases. 
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In principle one could check this by making a separate variational calculation for the 15 N 
states, but the sensitivity of the results on the shape of the s.p. wave functions is not 
large enough to do this reliably. In order to have an idea about the magnitude of this 
effect we calculated p-shell spectroscopic factors in a fixed-center mean-field model, using 
slightly different W.S. well shapes for the protons in 15 N and 16 0, and keeping the neutron 
wells identical. In this case it is easy to show that S p = Q^Q^ , with Q s (Q p ) the overlap 
between the (slightly) different proton s-orbitals (p-orbitals) in 16 and 15 N. If we fitted 
the W.S. parameters to the experimental charge densities available in the literature |I6 |, we 
find S p = 0.946, i.e. a substantial 5.4% reduction. However, a minimal change in the 16 
geometry (just reducing the W.S. radius to reproduce the rms radius of 15 N) gives only a 
1.2% reduction. We tend to regard the latter value as the more realistic. The rms radius 
of 15 N is better established than the whole charge density. The 15 N charge density in [16 



has the peculiarity of overshooting the 16 O charge density in the interior, and this feature is 
responsible for the more substantial effect. Also note that, lacking information, the neutron 
wells were taken identical in 15 N and 16 O. Relaxing this assumption would also contribute to 
the reduction of the proton p-shell spectroscopic factor. We estimate therefore the reduction 
of S P3/2 due to changes in the s.p. wave functions to be at most about 5% and probably less. 



B. Nonorthogonality and antisymmetrization corrections 

Here we look at the effects of nonorthogonality of the final scattering states ^(A)p i n 
Eq. ( |TUD and the antisymmetrization correction ( |TTD of the amplitude. With fixed-center wave 
functions these effects are known to be small under normal (e, e'p) kinematical conditions, 
but it is not inconceivable that cm. effects change this for intrinsic wave functions. Note 
e.g. that in the fixed-center frame the orthogonality corrections vanish if the final state is 
a SD with one of the bound hole states replaced with a continuum state of the same mean 
field, whereas this is no longer the case when ISD are used. 

In order to handle the nonorthogonality corrections on the amplitude, we consider the 
Lowdin transformation on the set of final scattering states ^(A)p W%- The orthonormal set 
of Lowdin transformed states can be expressed as 

l*(A)p> = / dAN- 1 ^ 2 ]^^), (28) 

in terms of the overlap matrix [N], 

[N] pp/ = (* {A)p ,\* (A )p) = S(p - p') - C(p, p'), (29) 

given by Eq.(|i"lj). It can be shown that this transformation does not change the correct 
asymptotic behavior (0) of the final scattering states. 

The amplitude can now be expanded in powers of the small correction C, and up to first 
order we get 

(*(A) P |p(q)|*(A)i) = (*(A) P |p(q)|*(A)i> + 5/ ^p / C(p,p , )(^(A)p'|p(q)|^(A)i). (30) 

Using Eqs. ( p^6] - p~7| ) and neglecting antisymmetrization corrections in the (small) second term 
of Eq . (|30D , the amplitude becomes 
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(*(A) P |p(q)|*(A)i) = Ap W iA(p,q) + A A .s.(p,q) + ^L6w(p,q), (31) 
with the first order Lowdin correction given by 

^L6w(p, q) = -- J dp'C(p, p / )A PWIA (p / , q). (32) 



In the harmonic oscillator model the magnitude of the different terms in Eq.(plf) can be 
easily estimated by concentrating on the Gaussian part of the momentum dependence. As 
an example, the exponent in the Gaussian part of the left-hand-side of Eq. (|5|) contains 

A A— 1 a i A 1 

5>i - Ra) 2 + J>* - Ra-i) 2 = — <-(rA - Ra) 2 + 2 £(r, - R^) 2 , (33) 
i=i i=i A i=i 



and it follows that 



1 A — 1 

ip(x) = exp[ — — - — r 2 ] x [Polynomial in r]. (34) 

AO J\ 



As a consequence, the Gaussian part of the momentum dependence of Apwia is 



-4pwia ~ ex P[ _ ^"23i( p ~ ~A~ q ^- ^ 



Similarly we find 



p{x,x') ~exp[-^(r 2 + r' 2 )], (36) 

C(p, p') ~ exp[-^^{(p' + -^p) 2 + (p + ^P') 2 }], (37) 

b 2 A , , A-l 2l 



/fa* A A 1 
dp'C(p, p) ex P[~YZ^T^ P ' A~ q 



~ exp[ 4 I3T {(P + I q)2 + ^/ q2}] - (38) 

The last result also holds for the antisymmetrization correction A A .s.(p, q)- 

The NIKHEF 16 0(e, e'p) experiment JTJ| was performed under quasi-elastic parallel kine- 
matics, with a roughly constant laboratory proton kinetic energy T PL « 90 MeV and missing 
momentum p m = Pl ~ Ql scanned in the region (-150 , 250) MeV/c. Under these kinemat- 
ical conditions the magnitude] 2 ] of the correction term ( |38| ) is very small compared to the 
magnitude (|35"D of the leading PWIA term in the amplitude, the ratio ranging from 4 x 10~ 6 
at the most negative p m , to 3 x 10~ 5 at p m =0, and growing to 6 x 10~ 3 at the most positive 
p m (corresponding to the smallest qi). To get a 10% correction at p m = 0, one would need 
protons ejected with only T PL m 20 MeV. We conclude that, under normal (e, e'p) kinemat- 
ics, the nonorthogonalityQ and antisymmetrization corrections are unimportant also for the 
translationally invariant wave functions considered here, just as in the fixed-center case. 



2 Note that in the laboratory frame the relative momentum p in Section II. B is given by p = 
PL ~ j[QL 

3 There is also the effect of nonorthogonality between the scattering states ^(A)p an d the initial 
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V. DISCUSSION 



The observed magnitude of the experimental cross sections of (e, e'p) reactions at small 
missing energies is not satisfactorily explained by present theoretical models. For 16 we 
pointed out that cm. correlations enhance the cross section leading to p-shell quasi-hole 
states by about 7%. As a result theoretical predictions [Q-Q that consider effects primarily 
due to short-range and tensor NN correlations give S P3/2 ~ 0.97 for 15 N. Including also 
low-energy configuration mixing in the target and residual nucleus || provides additional 
depletion and fragmentation of strength and lowers this to S P3/2 ~ 0.81, which is still far 
above the experimental value S Ps . 2 (6.32) = 0.53 ± 0.05. At present this discrepancy does 
not seem to be understood. 

More generally, the large difference between the experimental value of the p-shell quasi- 
hole strength and the present variational result indicates that an important ingredient is 
missing in the variational wave function. This is also signaled by recent calculations for 
A < 7 nuclei ||18|| , in which it is found that the quality of the VMC wave function (used as 
input in subsequent Green's function Monte Carlo calculations) deteriorates with increasing 
A, when compared to the final GFMC result. The possibility of a-cluster components in 
the surface part of the wave function should be looked at. 
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APPENDIX A: HARMONIC OSCILLATOR MODEL FOR THREE PARTICLES 

IN A S 2 P CONFIGURATION 

Consider one spin-down and two spin-up neutrons in a harmonic oscillator well. This is 
the simplest non-trivial illustration of the h.o. model. The ground-state configuration has a 
spin-up and a spin-down neutron in the s-shell and one spin-up neutron in the p-shell. Its 
intrinsic wave function is 

b 2 3 

*(3)(ri,r 2 ,r 3 ; spins) = exp[- — ^(r, - R 3 ) 2 ] {{z x - Z 3 )(\uud) - \udu)) 

j'=i 

+ {z 2 — Z 3 )(\duu) — \uud)) + (z 3 — Z 3 )(\udu) — \duu))} , (Al) 



state which is straightforward to include in this analysis. The correction term to the ampli- 

tude has a momentum dependence exp[— ^-(xrrP 2 + icr^ 2 )] an d is of the same order of magnitude 
as the other correction terms. 
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where the spin states are differentiated by u, d, the occupied p-orbital is chosen in the z- 
direction, and the h.o. length parameter is b. 
The three possible one-hole states are 

b 2 2 

*(2) p ( r i> r 2> spins) = exp[— - J2( r j ~ R 2 ) 2 ](M) - \du)), 

i=i 

b 2 2 

^(2)s( r i. r 2> s P ins ) = ex PhirII( r i ~ R 2) 2 ](0i - Z 2 )\ud) - (z 2 - Z 2 )\du)), 

A i=i 

*W( ri ,r 2 , spins) = expf-^fo ~ K 2 ) 2 }(zi - z 2 )\uu). (A2) 

In the notation of Eq.(|2^3D we have = «W = n («) = i ; ^(«) = 2 and A® = 1. 
The spectroscopic factors for the valence p^ and hole states are, according to Eq. (p2]) , 
'-'p^ = | and Sjfi = (|)° = 1. The spectroscopic factor for the hole state can be found 
from the sum rule (|3|), = 2 — § = |. 

The same results are easily found by applying Eqs.(§-^|) with the intrinsic wave functions 
T|- [A2|) . The normalizations of the wave functions are 



7T 3 3 5 / 2 



r 3r 

(*e»)l« r 0O> = L -Di-. (A3) 



b 8 



„3/2ol/2 

= !L f~. (AS) 



„3/2 9 3/2 

= (A6) 

The overlap functions are then given by 

fo 5/2o5/4 ft 2 

^(ra) = -v^3e*p[--r 3 2 ], (A7) 

A3/2ol/4 ft2 

V>f°(r.) = -^esp[--r 3 2 ], (AS) 

t3/2 2 3/4 ft 2 

^(rs) = ^ 7i exp[--r 3 2 ], (A9) 

and their normalization agrees with the values for the spectroscopic factors mentioned above. 

The spectroscopic sum rule (|23| ) in terms of the intrinsic one-hole states is a consequence 
of the special nature of a harmonic oscillator mean-field; in contrast to fixed-center systems 
it does not hold for general mean-field s.p. wave functions. As an example we can distort 
the h.o. mean field by taking different h.o. length parameters b and b' for the s and p orbitals 
in the present s 2 p model. The result for now becomes 

c( d ) _ 3 (^ 2 + Ux + 9)~ 7/2 [(x - l) 2 (x + 3) 2 + 16(x + l)(x 2 + Ux + 9)] 
3 ~ (x+l)[(2(2x + l))- 5 / 2 + 16(2(x + l)(x + 5))- 5 / 2 ] ' { } 
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with x = (j) 2 - For x = 1 we recover the result Sjfi = 1, i.e. in a pure h.o. mean field 
the spectroscopic strength of the spin-down neutron is fully contained in the intrinsic 
one-hole state. For x 7^ 1 we find < 1; the remainder of the strength is contained in 
more complicated configurations. 
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